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ABSTRACT 

We propose a phenomenological model for incompressible magnetohydro dy¬ 
namic turbulence. We argue that nonlinear-wave interaction weakens as the 
energy cascade proceeds to small scales, however, the anisotropy of fluctuations 
along the large-scale magnetic field increases, which makes turbulence strong at 
all scales. To explain the weakening of the interaction, we propose that small- 
scale fluctuations of the velocity and magnetic fields become increasingly dynam¬ 
ically aligned as their scale decreases, so that turbulent “eddies” become locally 
anisotropic in the plane perpendicular to the large-scale magnetic field. In the 
limit of weak anisotropy, that is, weak large-scale magnetic field, our model re¬ 
produces the Goldreich-Sridhar spectrum, while the limit of strong anisotropy, 
that is, strong large-scale magnetic field, corresponds to the Iroshnikov-Kraichnan 
scaling of the spectrum. This is in good agreement with recent numerical results. 

Subject headings: MHD-turbulcnce 


1. Introduction. 

Magnetohydrodynamic (MHD) turbulence is a state of a randomly stirred conducting 
fluid in the limit of very small fluid viscosity and resistivity; it plays an essential role in 
a variety of astrophysical systems, from planets and stars, to interstellar and intergalactic 
media (see, e.g., Biskamp 2003). Despite more than 35 years of analytical, numerical, and 
observational investigations, the spectrum of MHD turbulence remains a subject of contro¬ 
versy. The standard results and recent developments in the theory of MHD turbulence can 
be found in many excellent texts (see, e.g., Iroshnikov 1963; Kraichnan 1965; Goldreich & 
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Sridhar 1995, 1997; Ng & Bhattacharjee 1996; Clio & Vishniac 2000; Biskamp & Miiller 2000; 
Maron & Goldreich 2001; Clio, Lazarian & Vishniac 2002; Galtier, et al. 2000, 2002; Miiller, 
Biskamp & Grappin 2003; Biskamp 2003; Kulsrud 2005). In the present Letter we first an¬ 
alyze the approaches of Iroshnikov (1963) and Kraichnan (1965) and Goldreich & Sridhar 
(1995) and point out the discrepancies of these theories with numerical results. Then, we 
propose a new model for the turbulent MHD cascade, which is free of these discrepancies, 
and which is in good agreement with recent high-resolution numerical findings of Maron & 
Goldreich (2001), Miiller, Biskamp & Grappin (2003), and Padoan, et al. (2004). Our results 
are summarized in Conclusions. 


The Iroshnikov-Kraichnan energy cascade. The MHD equations describing evolution 
of the fluid velocity field, v(r, t), and fluctuations of the magnetic field, b(r,f), have an 


especially simple form when expressed in the Elsasser variables, z = v — b, and w = v + b: 

8 t z + (V A • V)z + (w • V)z = -VP, (1) 

d t w — (V A ■ V)w + (z • V)w = —VP, (2) 

where the pressure P is determined from the incompressibility condition, V • z = 0 or 


V • w = 0. In equations (1) and (2), V A = B 0 /v^rp is the Alfven velocity, p is the fluid 
density, Bo is the large-scale external magnetic field, and we have neglected small viscosity 
and resistivity. Note that for w = 0, any function z = f(r — V A f) is a solution of the 
system (eqs. 1,2); analogously, for z = 0, the solution is w = g(r + V A f), where the 
function g is arbitrary. 

Iroshnikov and Kraichnan used this fact to propose that the interacting Alfven-wave 
packets (or “eddies”) are those propagating in opposite directions along the large-scale 
magnetic-field lines (Iroshnikov 1963; Kraichnan 1965). From that, they deduced that the 
energy transfer time from wave packets of size A to smaller ones is increased compared with 
the simple dimensional estimate r(A) ~ X/Sv\. (We denote by Sv\ and 5b\ the velocity and 
magnetic-field fluctuations in the “eddy.” In the Alfven wave, Sv\ ~ 5b\.) Indeed, consider 
two wave packets of size A propagating in the opposite directions along a magnetic-field line 
with the Alfven velocities. Assuming that the “eddies” are decorrelated at the field-parallel 
scale A, one can estimate from equations (1) and (2) that during one collision the distortion 
of the “eddy” is A 5v\ ~ (5v\/ X)(X/Va). The distortions add up randomly, therefore, the 
“eddy” will be distorted relatively strongly after N ~ (5v\/A5v\) 2 ~ (Va/Sv \) 2 collisions. 
The energy transfer time is, therefore, 

tik( A) ~ NX/V a ~ X/6v x (y A /Sv x ). (3) 


It is important that in the Iroshnikov-Kraichnan interpretation, an “eddy” has to ex- 
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perience many uncorrelated interactions with oppositely moving “eddies” before its energy 
is transferred to a smaller scale. Moreover, turbulence becomes progressively weaker as the 
energy cascade proceeds toward smaller scales. The requirement of constant energy flux over 
scales = const immediately leads to the scaling of fluid fluctuations Sv\ oc A 1 / 4 , 

which resulted in the energy spectrum (the Iroshnikov-Kraichnan spectrum), 

E IK (k) = \8v k \ 2 k 2 oc k~ 3/2 . (4) 

Iroshnikov and Kraichnan did not consider anisotropy of the spectrum, so Eix(k) was as¬ 
sumed to be three-dimensional and isotropic. 


Anisotropy: The Goldreich-Sridhar energy cascade. Over the years, isotropy of the MffD 
spectrum in a strong external magnetic field seemed to contradict analytical and numerical 
findings, (see, e. g., Biskamp 2003). A theory of anisotropic MHD turbulence was proposed 
by Goldreich & Sridhar (1995). They argued that “eddies” are strongly anisotropic; they are 
elongated along the large-scale magnetic field lines. As a consequence, the time of relatively 
strong “eddy” distortion (the energy transfer time) is on the order of a crossing time required 
for two oppositely moving “eddies” to pass each other. 

Suppose that the “eddy” has a transverse (to the large-scale field) size A. Then, its field- 
parallel size l can be found from the so-called critical-balance condition, proposed by Gol¬ 
dreich & Sridhar (1995). This condition has two explanations that are equivalent in the 
Goldreich-Sridhar picture. First, the critical balance can be understood as a formal balance 
of the second and third terms in equations ( 1 ) and ( 2 ): Va/1 ~ 8b\/ A. Second, it follows 
from the geometric distortion of magnetic-field lines in the turbulent “eddy.” Indeed, the 
“eddy” displaces the lines in their perpendicular direction by a distance £ ~ 5b\l/VA, and 
this displacement equals the perpendicular “eddy” size A. 

The critical balance condition is the same for all scales, so, contrary to the Iroshnikov- 
Kraichnan picture, the turbulence strength does not change with the scale. The energy 
transfer time predicted in the Goldreich-Sridhar theory is 

t G s( A) ~ l/V 4 ~ A /8v x . (5) 

Assuming that the energy cascade is independent of the scale, 8v\/tgs{ A) = const, one 
obtains the scaling of fluid fluctuations, 5v\ oc A 1 / 3 . The corresponding energy spectrum is 

E G s(k±) = \8v k± \ 2 k± (x k^_ 5/3 . ( 6 ) 

This spectrum coincides with the Kolmogorov spectrum of incompressible non-magnetized 
fluid turbulence, as it should since the energy transfer time coincides with the “eddy” turn¬ 
over time, r(A) ~ A /5v\, in both approaches. The anisotropy of fluctuations is described by 
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the condition that follows from the critical balance, l oc A 2 / 3 . One can therefore write that 
the fluctuations scale with the field-parallel size of the “eddy” as Sv[ oc Z 1 / 2 . 

The Goldreich-Sridhar picture, however, does not fully agree with numerical simulations. 
As was recently discovered by M filler, Biskamp & Grappin (2003), the anisotropic spectrum 
depends on the strength of the external magnetic held. Denote 7 = /{p8v \), where 5vl is 

the velocity held at the outer scale of turbulence, L. ft was found that the field-perpendicular 
scaling of fluctuations changed from the Goldreich-Sridhar form to the Iroshnikov-Kraichnan 
form as the held increased from 7 < 1 to 7 > 1. A similar result for 7 3 > 1 was obtained 
earlier by Maron & Goldreich (2001). These intriguing numerical findings motivated our 
interest in the problem. In the next section, we propose a phenomenological model of MHD 
turbulence, which agrees well with available numerical results, for any strength of the external 
held. In the limiting case of a weak external held, our model reproduces the anisotropic 
spectrum of Goldreich & Sridhar. In the other limiting case of a very strong external held, 
anisotropy of the spectrum is stronger, however, the field-perpendicular spectrum formally 
coincides with the spectrum predicted in the Iroshnikov-Kraichnan model. 


2. A model for MHD turbulence. 

To begin with, we make a certain assumption about reduction of the nonlinear interac¬ 
tion, which is not present in either the Iroshnikov-Kraichnan or Goldreich-Sridhar picture. 
We postpone the justihcation of this assumption until the end of this section, when we ob¬ 
tain the corresponding solution for the turbulent spectra and compare it with numerical 
simulations. 


Reduction of nonlinear interaction. Let us assume that the nonlinear interaction of 
the counter-propagating fluctuations is reduced by a factor (<5u A /VA) a , where a is some 
undetermined exponent, 0 < a < 1. In other words, we assume that the nonlinearity in 
equations ( 1 ) and ( 2 ) is “depleted,” so that the interaction terms are of order 

(w • V)z ~ (z • V)w ~ [fv^lX){fv\!VA) a . (7) 

Thus, the fluid fluctuations at the transverse distance A become decorrelated on the time- 
scale 


M\)~(\/8v x )(V A /5v x ) a . ( 8 ) 

Their decorrelation length along the magnetic held line can be found from the causality 
principle. For 5v\ < Va, the perturbation cannot propagate along the held line faster than 
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Va; therefore, the correlation length along the held line is on the order of l ~ VAtjv(A). 
This condition is analogous to the critical balance condition of the previous section in that it 
satisfies the same balance between the linear and nonlinear interaction terms in equations (1) 
and (2). However, its geometric meaning is different, and will be discussed below. 

We see that the interaction strength (eq. 7) decreases for smaller scales; however, the 
field-parallel “eddy” size l is adjusted in such a way that the energy transfer to the smaller 
scales always takes one “eddy” crossing time. Note that contrary to the Iroshnikov-Kraichnan 
formalism, and similar to the Goldreich-Sridhar approach, in our picture turbulence is strong 
and fluctuations are highly anisotropic. Now, we require that the energy cascade be constant 
over scales 5 v\/tn(\) = const. We derive 8v\ oc A 1 ^ 3 " 1 ""), and the anisotropy condition reads 
l oc A 2 /( 3+ “). The corresponding energy spectrum is given by 

E(k ± ) = \8v k± \ 2 k ± oc kl {5+am3+a) . (9) 

One can formally define the corresponding longitudinal spectrum of fluctuations, from the 
condition E(k±)dk± = E(k\\)dk\\, with k± oc k^ +a ^ 2 . The answer is E(k\\) oc k ^ 2 ; note 
that it is independent of a. Therefore, the scaling of fluid fluctuations with respect to the 
field-parallel distance / is always Svi oc l 1 ^ 2 . 


Comparison with numerical simulations. Obviously, our result with a = 0 corresponds 
to the Goldreich-Sridhar scaling, while a = 1 produces the Iroshnikov-Kraichnan scaling. 
Simulations by Milllcr, Biskamp & Grappin (2003) for a range of large-scale external magnetic 
fields have shown that the scaling of the second-order structure function with respect to the 
field-perpendicular scale A changes from the Goldreich-Sridhar value in the case of a weak 
external field, 7 -C 1, to the Iroshnikov-Kraichnan value in the case of a strong field, 7 3> 1. 
At the same time, the scaling of the second-order structure function with respect to the 
field-parallel distance l does not change much and stays close to {(8vi) 2 ) oc l 0 ' 9 . This result is 
consistent with our prediction, (8vi) 2 oc l. The same scalings for the case of a strong external 
magnetic field, were earlier observed in simulations by Maron & Goldreich (2001). In these 
simulations, it was also found that the scaling of the energy-cascade time was te oc A 1 / 2 . 
This contradicts the Goldreich-Sridhar picture (see the discussion in § 6.1.3 of Maron & 
Goldreich 2001), but coincides with our formula Tat(A) oc A^ 1+a ^^ 3+Q \ for a = l. 

It is important to note that when the perpendicular Iroshnikov-Kraichnan scaling is 
formally reproduced in our model, the turbulent fluctuations are strong and anisotropic, 
with l oc A 1 / 2 . Moreover, the turbulence possessing the Iroshnikov-Kraichnan spectrum is 
more anisotropic than the turbulence in the Goldreich-Sridhar model, which is not surprising, 
since, as just discussed, the Iroshnikov-Kraichnan spectrum corresponds to a much stronger 
external magnetic field. 



6 


Scale-dependent dynamic alignment. We next explain a possible physical origin for the 
proposed reduction of the non-linear interaction (7). A hint toward the explanation can 
be obtained from geometric considerations. The displacement of magnetic held lines in the 
direction perpendicular to the large-scale magnetic held, produced by the wave packet with 
l oc \ 2 S 3+a l and 5b\ oc A 1 ^ 3 " 1 ""), is given by £ ex 5b\l ex A 3 /( 3+q A We thus obtain that 
the transverse displacement of magnetic held lines in the shear-wave packet is on the order 
of £ at distances A <C £. Therefore, this packet should be highly anisotropic in the plane 
perpendicular to the large-scale magnetic held. In this plane, it is elongated in the direction 
of the held huctuations. This may be consistent with the numerically supported picture that 
dissipative structures in MHD turbulence are micro-current sheets (Biskamp 1993, 2003; 
Biskamp & Miiller 2000; Maron & Goldreich 2001). For comparison, the Goldreich-Sridhar 
scaling, a = 0, would correspond to hlament-like dissipative structures. In Figure 1 we 
sketch the shapes of the “eddies” in the Goldreich-Sridhar picture and in our picture. In 
our model, the anisotropy of the “eddies” is weakest for a = 0, and therefore it is natural 
to identify this limit with zero large-scale magnetic held. In the other limit, a — 1, the 
anisotropy of the “eddies” is strongest, therefore, this case should be identified with a strong 
magnetic held, 7 > 1. 1 

Let us now explain how the envisioned anisotropy of huctuations in the held-perpendicular 
plane reduces the strength of the nonlinear-wave interaction. For this purpose, consider the 
third terms in the MHD equations (1,2), say, (w-V)z. Since both w and z are divergence-free 
and perpendicular to the large-scale magnetic held, this term is proportional to the angle 9\ 
between the directions of w\ and z\ (if this angle is small). But we just established that this 
angle is 9\ ~ A/£ oc A Q / ( ' 3+ "\ Quite remarkably, we have reproduced the reduction factor 
( 5v\/Va ) a in equation (7). This demonstrates that our initial assumption is self-consistent. 
Moreover, a detailed analysis of equations (1) and (2) suggests that the alignment 9\ oc 6v% 
is indeed consistent with the MHD dynamics, when 6\ is small. In this limit we can obtain 
the following approximate evolution equations, d t 5v\ rs_/ 8v\0\/\ and d t 9\ ~ 8v\9\/\. The 
alignment exponent a would be determined by numeric coefficients in these equations, which 
cannot be obtained from the dimensional analysis. 

As an important analogy, we note that decaying MHD turbulence approaches the so- 
called Alfvenic state, where either w or z is zero depending on the initial conditions (see 
e.g., Grappin, et al 1982; Grappin, Pouquet, & Leorat 1983). Based on our analysis, we 
propose that driven turbulence behaves in a similar manner, although b does not become 


■'■As is shown in (Grappin, et al 1982) with the aid of the EDQNM closure, strong large-scale magnetic 
field leads to the Alfvenic decorrelation time in aligned turbulence. In our model, this means that a cannot 
exceed 1. 
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exactly equal to ±v. Rather, magnetic fluctuations <R)a tend to align their direction, but not 
their magnitude, with that of the velocity fluctuations, ±5v A . As a result, the fluctuations 
5w\ and 5z\ are of the same order, and the directions of 5w\ and ±5 z\ are aligned within 
the angle 9\ ~ A/£ oc A Q ^ 3+a A The degree of the alignment increases progressively as the 
scale decreases. Such scale-dependent dynamic alignment (and the associated depletion of 
nonlinearity) can in principle be checked numerically, although the numerical analysis may 
be complicated by the rather slow dependence of the alignment angle 9\ on the scale. There 
is, however, a numerical indication that MHD turbulence indeed has a tendency to create 
correlated regions of polarized fluctuations (Maron & Goldreich 2001). 

On the non-universality of the turbulent spectrum. Our analysis points to an interest¬ 
ing possibility that MHD turbulence is non-universal in that it depends on the large-scale 
magnetic field. We may further speculate that, in principle, other large-scale conditions may 
affect the scaling properties of turbulence. For example, the dynamic alignment may be 
sensitive to the level of cross-helicity fluctuations; an analogous result is known for the case 
of decaying turbulence, (see, e.g., Grappin, et al 1982; Grappin, Pouquet, & Leorat 1983). 
An alternative possibility is that the spectrum is universal and has the Iroshnikov-Kraichnan 
scaling but, for 7 -C 1, the dynamic alignment with a = 1 is established rather slowly as the 
scale decreases, and the resolution of numerical simulations is not large enough to reach the 
universal regime. 


3. Conclusions. 

To conclude, we summarize our main results: 

1. Our consideration is motivated by the recent numerical observations (Maron & 
Goldreich 2001; Muller, Biskamp & Grappin 2003) that incompressible MHD turbulence is 
not completely described by either the Iroshnikov-Kraichnan or the Goldreich-Sridhar model. 
The scaling of velocity fluctuations was found in these papers to depend on the strength of 
the large-scale external magnetic field: the Iroshnikov-Kraichnan scaling appeared in the 
limit of strong external magnetic field, while the Goldreich-Sridhar scaling appeared in the 
limit of weak field. 

2. To explain these numerical findings, we propose that turbulent fluctuations become 
increasingly dynamically aligned as the energy cascade proceeds to smaller scales. Velocity 
fluctuations Sv\ tend to align their direction with that of magnetic field fluctuations, ±<5bA, 
and the smaller the scale A, the stronger the alignment. The dynamic alignment leads to 
reduction of the nonlinear-wave interaction (so-called depletion of nonlinearity). 



3. As a result of point 2 , fluctuating “eddies” are three-dimensionally anisotropic. The 
“eddies,” whose smallest scale is A, have the scale £ oc A 3 ^ 3 " 1 " 0 -* in the direction of the shear 
(the direction of the magnetic field line distortion), and the scale l oc A 2// ^ 3+ “^ in the direction 
of the large-scale magnetic field, as is sketched in Figure 1. The scaling and anisotropy of 
fluctuations are described by a single parameter a, which depends on the strength of the 
external magnetic field, and which is determined in our work from comparison with numerical 
simulations by Maron & Goldreich (2001) and Miiller, Biskamp & Grappin (2003). 

4. The energy distribution is given by 

E(k ± ) OC £;-( 5 +«)/(3+«). 

This coincides with the Goldreich-Sridhar spectrum in the limit of weak anisotropy (a = 0), 
and with the Iroshnikov-Kraichnan spectrum in the limit of strong anisotropy (a = 1). 
As the external magnetic field increases from 7 < 1 to 7 > 1 , the anisotropy of turbulent 
fluctuations increases, and the parameter a changes from a = 0 to a = 1. According to point 
3, the corresponding scalings of the fluctuations with respect to their field-perpendicular and 
field-parallel dimensions are Sv\ oc A 1 ^ 3 " 1 ""), Jty oc £ 1//3 , and Svi oc l 1 ' 2 . 

5. The smallest-scale “eddies” in our turbulent cascade (A —» 0 for a ^ 0) have a sheet 
like morphology, in agreement with micro current-sheet dissipative structures, numerically 
observed in MHD turbulence (Biskamp & Miiller 2000). In the case of zero external mag¬ 
netic field, a = 0 , the dissipative structures are filaments, which also agrees with numerical 
simulations (Padoan, et al. 2004). 

6 . Previous attempts to explain the numerically observed Iroshnikov-Kraichnan scaling 
(and the associated cascade-time increase (3)), in anisotropic sub-Alfvenic MHD turbulence 
(7 3> 1) essentially invoked intermittency effects, (see, e.g., § 6.6.4 in Maron & Goldreich 
2001; Muller, Biskamp & Grappin 2003). The intermittency effects are essential for ex¬ 
plaining higher order statistics of MHD turbulence, that is, scaling of higher order structure 
functions of w and z. However, they usually provide only small corrections to the scaling 
of the second-order structure functions, and of the energy spectra. Such effects are not ad¬ 
dressed in the present Letter; our derivation is based on the idea of scale-dependent dynamic 
alignment, which does not require intermittency. 
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Fig. 1.— A sketch of the shapes of turbulent fluctuations, “eddies,” in the Goldreich-Sridhar 
picture (left), and in our picture (right). The large-scale magnetic held is denoted by B. As 
the turbulent cascade proceeds toward the smallest, dissipative scales, A —> 0, the Goldreich- 
Sridhar “eddy” assumes the shape of a filament, while the “eddy” predicted in our model 
turns into a current sheet. 
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